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Abstract 

We simulate lattice models in the three-dimensional XY universality class 
in the low and the high temperature phase. This allows us to compute a 
number of universal amplitude ratios with unprecedented precision: i?x = 
0.411(2), R B = 2.83(1), = 0.3562(10) and RJ = 0.850(5). These results 
can be compared with those obtained from other theoretical methods, such 
as field theoretic methods or the high temperature series expansion and also 
with experimental results for the A-transition of 4 He. In addition to the 
XY model, we study the three-dimensional two-component 4 model on the 
simple cubic lattice. The parameter of the 4 model is chosen such that 
leading corrections to scaling are small. 
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1 Introduction 



In the neighbourhood of a second order phase transition various quantities diverge, 
following power laws. E.g. in a magnetic system, the magnetic susceptibility be- 
haves as 



where 7 is the critical exponent of the magnetic susceptibility, C + ,C_ are the 
amplitudes in the high and low temperature phase, respectively, and t = (T — T c )/T c 
is the reduced temperature. Critical exponents like 7 assume universal values; i.e. 
they assume exactly the same value for all systems within a given universality 
class. Following Wilson (see e.g.p]), such a universality class is characterized by 
the dimension of the system, the range of the interaction and the symmetry of the 
order parameter. In addition to the critical exponents, amplitude ratios like C + / C- 
are universal, while the value of C + or C_ depends on the microscopic details of 
the model. For a review on amplitude ratios see ref. [2]. 

In the present work we compute, using data obtained from Monte Carlo simula- 
tions of lattice models, the numerical values of four such universal amplitude ratios 
for the universality class of the three dimensional XY model. The A-transition of 
4 He is supposed to share this universality class. At temperatures below the tran- 
sition, 4 He becomes superfluid. The A-transition owes its name to the fact that 
the specific heat plotted as a function of temperature resembles the Greek letter 
A. The order parameter of the A-transition in 4 He is the phase of a wave func- 
tion. Therefore it should share the XY universality class, which is characterized 
by the 0(2), or equivalently U(l), symmetry of the order parameter. The study of 
the A-transition provides by far the most precise experimental results for universal 
quantities like critical exponents and amplitude ratios. Thus this transition gives us 
a unique opportunity to test the ideas of the renormalization group and to bench- 
mark theoretical methods. For a review and an outlook on future experiments in 
space-stations B see ref. [5] . 

The present work is the completion of ref. [6] , where we had computed the ratio 
A+/A- of the amplitudes of the specific heat. Our results can be compared with 
those obtained by using other theoretical methods and, which may be even more 
important, experimental results obtained for the A-transition of 4 He. 

The outline of our paper is the following: First we define the models and the 
observables that are measured. We briefly discuss the Monte Carlo algorithm that 
has been used for the simulation. Since the continuous 0(2)-symmetry is sponta- 
neously broken in the low temperature phase, there is a Goldstone boson. As a 
consequence, the thermodynamic limit is approached with corrections that decay 
as inverse powers of the linear lattice size. Therefore, extracting estimates for the 
thermodynamic limit from Monte Carlo data requires special care. To this end, we 

*The condition of micro-gravity avoids the broadening of the transition due to the gravitational 
field and hence allows to access reduced temperatures down to 5 x 10~ 10 [Sid]. 
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summarize the relevant results of chiral perturbation theory as discussed in refs. 
[3 [8j [9J . Then we present our numerical estimates of various observables in the low 
and high temperature phases. Next we compute the amplitude ratios from these 
observables. Finally we compare our results with those of other theoretical methods 
and experiments. 



2 The models 

We study the 4 model on a simple cubic lattice with periodic boundary conditions 
in each of the directions. The classical Hamiltonian is given by 

Hp = -0 £ & • + £ [<S + A(<g. - l) 2 ] , (2) 

<x,y> x 

where the field variable <fi x is a vector with two real components. < x,y > denotes a 
pair of nearest neighbour sites. The sites of the lattice are labelled by x = (xi,X2, £3) 
with Xi G {1, 2, Li}. Throughout we consider lattices with L\ = L 2 = L 3 = 
L. Note that in our convention, following ref. [ID], the inverse temperature /3 is 
absorbed into the Hamiltonian but does not multiply its second term. The partition 
function is given by 



n 



exp(-iy^) . (3) 



The reduced free energy density is then given by / = — -p log(Z). 

In the limit A — > 00 the classical XY (plane rotator) model is recovered. It is 
defined by the classical Hamiltonian 

H X y = -P Yl g * g y ' ( 4 ) 

<x,y> 

where s x is a unit-vector with two real components. For A = one gets the exactly 
solvable Gaussian model. For < A < 00 the 4 model undergoes a second order 
phase transition in the XY universality class; see e.g. ref. [XI] . 
Power laws like eq. Q are subject to corrections; see e.g. ref. 



x=c ± i*r 7 (i+M e + -) (5) 

with 9 pa 0.5 for the three-dimensional XY universality class. Such corrections 
complicate the determination of universal quantities from Monte Carlo simulations 
or high temperature series expansions of lattice models. Correction to scaling am- 
plitudes like b± depend on the parameter A of the 4 model. It is a rather old 
idea [12] to search for the value of such a parameter, where the leading correction 
to scaling amplitude vanishes. Note that the renormalization group predicts that 
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the zero of the leading correction amplitude is the same for all quantities. It has 
been demonstrated numerically [T3l [H] that such an improved value A* indeed ex- 
ists. The most accurate determination of the improved value is A* = 2.15(5) |10j . 
Previous estimates are A* = 2.07(5) in ref. [2] and A* = 2.10(6) in ref. [13J. 

We performed simulations at A = 2.1 and A = 2.2, where leading corrections 
to scaling are small. Following ref. [10] leading corrections to scaling at A = 2.1 
and A = 2.2 should be at least by a factor of 20 smaller than in the XY model. 
Comparing results obtained with A = 2.1 and A = 2.2 allows us to estimate the 
effect of these small corrections. In ref. [JU] the estimates (3 C = 0.5091503(3) [3] 
and j3 c = 0.5083355(3) [4] for the inverse critical temperatures at A = 2.1 and 
A = 2.2, respectively, were obtained. The number quoted in () gives the statistical 
error, while the number given in [] gives the systematic one. In the following 
analysis of our data we will just add up these two numbers. In order to clearly see 
effects due to leading corrections to scaling, we also have simulated the XY model. 
Recent estimates for the inverse of the critical temperature are f3 c = 0.454165(4), 
0.454167(4), 0.4541659(10) and 0.4541652(11) in refs. [151 M, EH HQ] , respectively. 
In the analysis of our numerical data we shall assume (3 C = 0.0.4541655(10), which 
is roughly the average of the estimate of refs. [17] and [TO] . 

2.1 Simulation algorithm 

We have simulated the XY model using the single cluster algorithm [T5] . The 4> 4 
model was simulated using a hybrid of the single cluster algorithm and a local 
Metropolis algorithm. The cluster algorithm allows only to change the angle of the 
field 0. In order to to change the modulus |0|, Metropolis updates are performed. 
Such a hybrid approach was originally proposed by Tamayo and Brower [19] for the 
one component 4 model. The generalization to iV components is discussed in ref. 
|13j . For details of our implementation, in particular of the local updates see ref. 

ma- 

3 Observables 

In the following we give the precise definitions of the quantities that we have mea- 
sured in our Monte Carlo simulations. To this end we use the notation of the 4 
model. The definitions for the XY model can be obtained by replacing <fi by s in 
these definitions. 

3.1 The energy and the specific heat 

In order to study universal quantities it is not crucial how the transition line in 
the j3-\ plane is crossed, as long as this path is smooth and not tangent to the 
transition line. Here, following computational convenience, we vary (3 at fixed A. 
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Correspondingly we define the energy density as the derivative of the reduced free 
energy with respect to (5. Furthermore we multiply by —1 to get positive numbers: 

tf = ^s( ( 6 ) 

\<x,y> I 

We then define the specific heat as the derivative of the energy density with respect 
to (3: 



L 3 







« 


£ <frx ■ <fry 




_<x,y> 




(7) 



3.2 The magnetisation, the magnetic susceptibility and the 
correlation length 

The magnetisation of a configuration is given by 



m 



V ^ 



The second moment correlation length is defined by 

^2nd — 



X/F-I 



4 sin(7r/L) 2 ' 

where the magnetic susceptibility in the high temperature phase is given by 



(9) 



(10) 



where we assume that the magnetisation vanishes in the high temperature phase. 
The Fourier transform of the correlation function at the lowest non-zero momentum 



is 



F k = — 
2V 



E 



exp l i 



In order to reduce the statistical error we have averaged F over all three directions 
k = 1,2,3. 

3.3 The helicity modulus T 

The helicity modulus T gives the reaction of the system under a torsion. To define 
the helicity modulus we introduce rotated boundary conditions in one direction: 
For Xi = Li and y x — 1 we replace the term <fr x <fr y in the Hamiltonian by 

0, • Rj y = 0« (cos(a)0« + sm(c04 2) ) + (fr ( ? (- sin(a)0« + cos(a)0f) (12) 
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The helicity modulus is then given by 



T 



L x <9 2 logZ(a) 



LoL 



2-^3 



do? 



(13) 



Q = 



Note that we have skipped for convenience a factor l/{ksT) compared with the 
standard definition. The helicity modulus can be directly evaluated in the Monte 
Carlo simulation. Following eq. (3) of ref. [20J one gets 



T 







L1L2L3 



-(1,0,0) 



P 2 



L1L2L3 



(,(2) 

x+(l,0,0) 



(1) 



(14) 

To arrive at this expression, the torsion by a is spread over the lattice ; i.e. using a 
torsion by a/L\ at any x\. In order to reduce the statistical error we have measured 
the helicity modulus also in X2 and ^-direction. 

The helicity modulus is of particular interest, since it plays a central role in the 
effective description of the behaviour in the low temperature phase, as we shall see 
below. Furthermore it can be accessed experimentally in the superfluid phase of 
4 He. In the literature, the inverse of the helicity modulus, as defined here, is called 
transverse correlation length. It is given by [21] 



1 _c 



m\k b T 
h 2 p s 



(15) 



where p s is the superfluid density and 7714 the mass of a 4 He atom. The superfluid 
density can be obtained from measurements of the second sound and the specific 
heat. 



4 The Goldstone mode and finite size effects 

Before we discuss the numerical results in the low temperature phase, let us summa- 
rize the relevant results from chiral perturbation theory. These results are derived 
in ref. [3, [8] for 0(iV)-invariant nonlinear o models in three dimensions. Let us 
briefly recall the results for our case N = 2, which is simpler than the general case. 
Furthermore, we consider only the case of a vanishing external field. 

The basic assumption of chiral perturbation theory is that in the broken phase, 
on scales much larger than the correlation length, only fluctuations perpendicular 
to the overall magnetisation remain as degrees of freedom. Their fluctuations can 
be described by an effective Hamiltonian H e ff[ip], where the field ip gives the angle 
of the fluctuation. Due to the 0(2) invariance of the underlying microscopic model, 
the effective model has to be invariant under global shifts of the field. I.e. it can 
only depend on derivatives of the field ip. Also by symmetry, it can only depend on 
even powers of odd derivatives. 



5 



In the leading, Gaussian approximation, the effective Hamiltonian, using lattice 
notation, is given by 

# e //M = ^£(^-Vv) 2 - (16) 

<xy> 

In the limit (3 — > oo, the effective Hamiltonian provides an approximation of the XY 
model in a microscopic sense with {3 e ff = (3. However, in the neighbourhood of the 
phase transition, it serves (only) as an effective model. This means in particular 
that the relation between (3 e ff and (3 is a priori not known. Hence f3 e ff has to 
be determined from observables that characterize the behaviour of the system at 
large scales. To this end, the helicity modulus is particularly useful. Plugging in the 
effective Hamiltonian H e ff[ip} into definition (1131) one obtains (3 e ff = T. Corrections 
to this relation are due to higher order terms in the effective Hamiltonian. One gets 



T(L) = T(oo) + c T L~ 



(17) 



for the L-dependence of the helicity modulus. A similar relation holds for the energy 
density: 

E(L) = E(oo) + c E L~ 3 . (18) 

In the case of the two-point function, a non-trivial effect arises from the fact 
that the spin is a non-linear function of the angle. This has to be taken into account 
in the case of the expectation value of the square magnetisation. From eq. (2.18) 
of ref. [9] we read 

-pi + 2p 2 a 2 



where 



(m 2 ) 



Pi = l + \l3 1 a + ^{(3 2 l -2l3 2 )a 2 



P'2 



with fix = 0.225785, (3 2 = 0.010608 and 



a 



1 

TL 



(19) 
(20) 

(21) 



These results are derived from more general expressions given in section 10 of ref. 
[7]. Corrections to eq. (fT9l) are proportional to a 3 . 



5 Simulations in the low temperature phase 

We have simulated the 4 model at A = 2.1 and A = 2.2 and the XY model for 
various values of (3 in the low temperature phase. Most of these simulations were 
already carried out in the context of ref. [6], where only the data for the energy 
density were analysed. For each value of (3 and each model we have simulated a 
number of lattice sizes to extrapolate, following the predictions of chiral pertur- 
bation theory, to the thermodynamic limit. In the case of the 4 model, for our 
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Table 1: Results for various lattice sizes L for the 2-component 4 model at A = 2.1 
and (3 = 0.515. We give numerical estimates for the helicity modulus T, the energy 
density E and the square of the magnetisation < m 2 >. S 2 is defined by eq. Alty). 
For a detailed discussion see the text. 



L 


T 


E 


(m 2 ) 


S 2 /2 


32 


0.05097(14) 


1.01763(8) 


0.09038(8) 


0.07822(7) [2] 


48 


0.04989(14) 


1.01655(5) 


0.08504(6) 


0.07725(5) [2] 


64 


0.04968(14) 


1.01620(3) 


0.08249(5) 


0.07677(5) [1] 


96 


0.04938(15) 


1.01606(2) 


0.08036(6) 


0.07661(6)[1] 


128 


0.04938(14) 


1.01603(1) 


0.07938(4) 


0.07657(4) [1] 



smallest values of j3, we have simulated lattices up to the linear size L = 288. Typ- 
ically we performed 10 5 to 2 x 10 5 measurements for each value of A, f3 and L. For 
each of these measurements one Metropolis sweep and a few single cluster updates 
were performed. In total the simulations in the low temperature phase took about 
1.5 years of CPU time on a 2 GHz Opteron CPU. 

As a typical example, we give in tabled] the results for A = 2.1 and (3 = 0.515. 
First we have fitted the helicity modulus T and the energy density E with the 
ansatze (I17|ll8j) . respectively. Using the data for all lattice sizes given in table [lj we 
get T(oo) = 0.049394(8) and E(oo) = 1.016005(8). These values are close to those 
obtained for the largest lattice size L = 128 that we have simulated for A = 2.1 and 
f3 = 0.515. Therefore we regard the extrapolation as save. In table [2] we summarize 
our results for T(oo) and E(oo) for all values of (3 that we have simulated in the low 
temperature phase of the 4 model at A = 2.1. Analogous results for A = 2.2 and the 
XY model can be found in tables [3] and HI respectively. In addition to the statistical 
error, we give in parenthesis the differences T(oo) — T(L max ) and E(oo) — E{L max ) 
as a rough estimate of the systematic error. L max is the largest lattice size that is 
available at a given value of (3. In order to compute the thermodynamic limit of 
the magnetisation we used the results of chiral perturbation theory summarized in 
section HI As input, we have taken our result for the thermodynamic limit of the 
helicity modulus. In the last column of table [U we give our results for S/2 for each 
value of L using eq. (fl9l) . In order to get an estimate for the thermodynamic limit, 
we have fitted these results with the ansatz E(L) = £(oo) + C5L 3 . Results obtained 
this way, for all values of (3 that have been simulated, are given in the last column 
of table El [3] and H] for A = 2.1, A = 2.2 and the XY model, respectively. 

Our results for the helicity modulus for the XY model can be compared with 
those of refs. 0, [TH]. We have taken the values of F from table 6 of ref. [TB]. The 
relation with the helicity modulus, as defined here, is T = F 2 . In the case of ref. [H] 
we have taken the same quantity from their table V. In figure [1] we have plotted the 
inverse of the helicity modulus as a function of (j3 c — j3)~ v using the numerical value 
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Table 2: Results for the low temperature phase of the 2-component <f) A model at 
A = 2.1. (3 is the inverse temperature, T the helicity modulus, E is the energy and 
T? the square of the magnetisation. For the dehnition of these quantities and a 
detailed discussion see the text. 



p 


T 


E 


£ 2 


0.51 


0.01330(15 


)[+4] 


0.931674(14) [-16] 


0.04020(5)[-7] 


0.5105 


0.01818(18 


)[-31] 


0.941232(16) [-13] 


0.05545(5)[-5] 


0.511 


0.02265(17 


)[-25] 


0.950382(11) [-9] 


0.06905(5) [-4] 


0.512 


0.03036(10 


)[-17] 


0.967852(11) [-24] 


0.09322(5)[-6] 


0.513 


0.03717(11 


)[-29] 


0.984474(10) [-29] 


0.11480(7) [0] 


0.515 


0.04939(8)| 


+ 1] 


1.016005(8) [-27] 


0.15305(7) [-12] 


0.52 


0.07456(7) 


-13] 


1.088210(8) [-17] 


0.23338(8)[+4] 


0.525 


0.09628(7) 


-3] 


1.15416(1) [-2] 


0.30099(9)[-14] 


0.53 


0.11565(7) 


-15] 


1.21557(2) [-7] 


0.36133(11) [-16] 


0.535 


0.13349(6) 


+6] 


1.27342(2) [-5] 


0.41614(11) [+2] 


0.54 


0.15035(6) 


-7] 


1.32830(2) [-4] 


0.46694(14) [+14] 


0.55 


0.18145(6) 


-8] 


1.43072(3) [-7] 


0.55824(17) [-33] 


0.58 


0.26340(5) 


-18] 


1.69533(6) [-35] 


0.78325(22) [-52] 



Table 3: Same quantities as in tabled for A = 2.2. 



p 


T 


E 


S 2 


0.5095 
0.51 
0.511 
0.512 


0.01670(16) [+2] 
0.02131(17)[0] 
0.02877(16) [-20] 
0.03602(24) [+17] 


0.937817(19) [-18] 
0.947000(19) [+5] 
0.964370(17) [-14] 
0.980920(19) [-14] 


0.04986(5) [-6] 
0.06386(5) [-6] 
0.08840(4) [-3] 
0.11022(5) [-4] 
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Table 4: Same quantities as in table [2| for the standard XY model. 



p 


T 


E 


E 2 


0.456 


0.02043(16) [-11] 


1.01820(1) 


[-2] 


0.06674(6) [-8] 


0.458 


0.03331(16)[-7] 


1.04568(2) 


[-3] 


0.11036(6) [-9] 


0.46 


0.04399(9) [-22] 


1.07095(1) 


[-6] 


0.14648(4) [-16] 


0.462 


0.05334(13)[-18] 


1.09473(2) 


[-15] 


0.17851(12)[-31] 


0.465 


0.06618(12) [+4] 


1.12832(2) 


[-9] 


0.22097(11)[-18] 


0.47 


0.08460(7) [+5] 


1.17991(1) 


[-9] 


0.28224(12) [-24] 


0.48 


0.11591(7) [-25] 


1.27155(5) 


[-27] 


0.38405(15) [-20] 


0.50 


0.16644(8) [-8] 


1.42298(3) 


[-21] 


0.53885(16)[-31] 


0.52 


0.20886(5)[-ll] 


1.54594(2) 


[-17] 


0.65719(18)[-2] 


0.525 


0.21859(6) [-2] 


1.57336(3) 


[-10] 


0.68335(4) [-12] 


0.55 


0.26344(6) [-14] 


1.69440(4) 


[-24] 


0.79552(39) [-10] 



v = 0.6717 for the exponent. The results of ref. [9] are in reasonable agreement with 
ours. On the other hand, those of ref. [16] clearly deviate. In particular for small 
reduced temperatures, the value of the helicity modulus seems to be underestimated 
in ref. [16] . Both ref. [5] and [IB] did not directly measure the helicity modulus 
using eq. ffl4|) but instead extracted it from the L-dependence of the magnetisation. 
In the case of ref. [16], in addition, the simulations were done for a finite external 
field, involving the extrapolation to vanishing external field. 

6 Simulations in the high temperature phase 

In the high temperature phase, we expect that the observables converge exponen- 
tially fast towards the thermodynamic limit. Therefore we have used lattice sizes 
that are large enough to ensure that the systematic error due to the finite lattice size 
is by far smaller than the statistical error of the observables that we have measured. 
We have checked that this is fulfilled for L ^ lO^nd- in fact for most of our the 
simulations we have chosen L > 12^2nd- Most of the simulations reported here were 
already performed in connection with refs. [TDJ [6]. The magnetic susceptibility x 
and F were measured using so called improved estimators; see e.g. ref. [22] • Note 
that these improved estimators are, in contrast to the naive ones, self- averaging. 
Hence, for large L/£ their statistical error is much smaller than that of the naive 
estimators. For the largest values of (3 we have performed 2 x 10 5 measurements. 
For each measurement we performed one Metropolis sweep and n c single cluster 
updates. The number n c of single cluster updates was chosen such that the average 
size of a cluster times n c is roughly equal to the number of sites of the lattice. The 
simulations for the largest lattice size L = 350 took about one month on a 2 GHz 
Opteron CPU for each value of A. Our results for the high temperature phase of 
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Figure 1: The inverse of the helicity modulus T is plotted as a function of (f3 c —(3)~ v 
using the numerical value v = 0.6717 for the exponent. Our results are given by 
the black circles; the data taken from ref. [9] are represented by red squares and 
the results of ref. [Tj6] are given by blue diamonds. The dotted lines should only 
guide the eye. 
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the 4 model at A = 2.1 and A = 2.2 and the XY model are summarized in tables 
[5j M and El respectively. 



7 Universal amplitude ratios 

In this section, we first define the amplitude ratios that we consider. Then we 
compute their values from the Monte Carlo data discussed in the previous sections. 
Finally we compare our results with those obtained by using other theoretical meth- 
ods and experimental results obtained for the A-transition of 4 He. 



7.1 Definition of the amplitude ratios 

Here we define the amplitude ratios that we have studied. For a more comprehen- 
sive list see e.g. page 18 of ref. [2] or table 2 of ref. [23J. The first amplitude 
combination relates the correlation length in the high temperature phase with the 
helicity modulus which is defined in the low temperature phase: 

Ry = /Wti+To , (22) 
where the amplitudes are defined by 

&nd=^f2ndA-tr U , T~T r , (t>0) (23) 

where v is the critical exponent of the correlation length. In this section we use, 
for computational convenience, 

t = P-P c (24) 

as definition of the reduced temperature. It has been shown in ref. [21] that 
the helicity modulus behaves as an inverse correlation length. The exponential 
correlation length £ exp , which describes the asymptotic decay of the correlation 
function, differs only slightly from £ 2 nd which has been used here. Following ref. 



Next we consider 



lim|^ = 1.000204(3) , (t < 0) . (25) 

Rb = j^j^ , (26) 

where the amplitudes of the magnetic susceptibility and the correlation length in 
the high temperature phase are defined by eqs. fllf23p . respectively. The amplitude 
of the magnetisation in the low temperature phase is defined by 

£/2~ (Bt p ) 2 , (t>0) . (27) 

As the third combination of amplitudes we consider 

R+ = M+) 1 ' 3 f 2ndt+ , (28) 
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Table 5: Results for the high temperature phase of the 2-component 4 model at 
A = 2.1. f3 is the inverse temperature, ^2nd the second moment correlation length, x 
the magnetic susceptibility, E the energy and C the specific heat. For the definition 
of these quantities and detailed discussion see the text. 
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0.40 


1 07551(2) 


5 7056(2) 


^1 268Q 


'10) 


2 091(3) 


0.41 


1 15844(2) 


6 4486(2) 


5S41 64 
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0.42 


1 25591(3) 
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8 6013(4) 
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'13) 


2 500 (4) 
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1 51582(4) 
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2 692(41 


45 
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1 697381' 41 


1 2 5045(5) 


635073 
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2 902(71 
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1 80821(4) 


14 0117(5) 
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3 050(81 


46 
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1 93718(5) 


15 8770(8) 
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'Si 


3 199(8) 


46^ 

W . TlW<J 


2 08960(6) 


1 8 2377(8) 


n fi89n9n 

U.UoZUZU 


o ) 


3 377f7) 


n 47 


2 27341 (8) 


21 3070(13) 


699359 


'9) 




0.475 


2.50035(11) 


25.4333(18) 


0.717754 


;q) 


3.813(8) 


0.48 


2.79047(12) 


31.2425(22) 


0.737430 


;q) 


4.081(10) 


0.485 


3.17678(14) 


39.909(3) 


0.758610 


!q) 


4.385(11) 


0.49 


3.72370(19) 


53.998(5) 


0.781682 


y] 


4.87(2) 


0.493 


4.1825(2) 


67.451(6) 


0.796679 


;e) 


5.09(3) 


0.495 


4.5763(2) 


80.177(6) 


0.807259 


|6) 


5.42(2) 


0.50 


6.1498(5) 


141.899(17) 


0.836367 


is) 


6.26(2) 


0.503 


8.0424(8) 


238.88(5) 


0.856373 


y] 


7.14(5) 


0.505 


10.4822(16) 


400.38(11) 


0.871351 


;§) 


7.95(3) 


0.506 


12.6264(16) 


575.74(14) 


0.879580 


;e) 


8.52(5) 


0.507 


16.318(4) 


950.7(4) 


0.888476 


7) 


9.43(8) 


0.5075 


19.498(6) 


1347.1(8) 


0.893283 


;q) 


9.83(7) 


0.508 


24.845(8) 


2164.6(1.4) 


0.898418 


!e) 


10.81(8) 


0.5083 


30.453(10) 


3225.0(2.0) 


0.901727 


I 4 ) 


11.32(9) 



Table 6: Same quantities as in tabled for A = 2.2. 



p 




2X 


E 


C 


0.501 


7.1723(5) 


191.372(23) 


0.849150(4) 


6.681(10) 


0.5035 


9.5018(9) 


330.80(6) 


0.866864(5) 


7.545(15) 


0.5055 


13.6104(21) 


667.11(20) 


0.882972(5) 


8.62(4) 


0.5067 


19.710(5) 


1376.2(6) 


0.894014(6) 


9.87(6) 


0.50748 


30.475(10) 


3231.4(2.0) 


0.902171(4) 


11.41(9) 



12 



Table 7: Same quantities as in table® for the standard XY model. 



p 


^2nd 


2X 


E 


C 


0.4 


1.87550(7) 


16.8465(10) 


0.740797(15) 


3.139(4) 


0.41 


2.18009(8) 


22.0757(13) 


0.773412(11) 


3.388(5) 


0.42 


2.62528(12) 


31.019(2) 


0.809006(12) 


3.736(6) 


0.425 


2.93916(12) 


38.242(2) 


0.828195(8) 


3.957(6) 


0.43 


3.35803(15) 


49.064(4) 


0.848624(8) 


4.231(7) 


0.435 


3.9514(2) 


66.702(5) 


0.870561(8) 


4.561(6) 


0.44 


4.8769(2) 


99.564(6) 


0.894452(5) 


5.014(8) 


0.441 


5.1305(2) 


109.704(8) 


0.899537(6) 


5.139(7) 


0.442 


5.4184(2) 


121.801(8) 


0.904735(5) 


5.242(8) 


0.443 


5.7489(3) 


136.464(12) 


0.910056(6) 


5.390(7) 


0.444 


6.1327(3) 


154.524(14) 


0.915514(6) 


5.542(7) 


0.445 


6.5851(4) 


177.23(2) 


0.921133(6) 


5.704(8) 


0.446 


7.1286(5) 


206.53(2) 


0.926937(6) 


5.874(10) 


0.447 


7.7960(4) 


245.53(2) 


0.932911(4) 


6.104(9) 


0.448 


8.6402(4) 


299.62(2) 


0.939131(3) 


6.342(9) 


0.449 


9.7488(6) 


378.69(4) 


0.945605(4) 


6.644(10) 


0.45 


11.2871(7) 


503.40(5) 


0.952402(3) 


6.959(10) 



where the amplitude of the specific heat Ch is given by 

C h ~ A±\t\~ a + b . (29) 

The analytic background b has to be taken into account, since the exponent of the 
specific heat a is negative for the three-dimensional XY universality class. The 
fourth combination of amplitudes that we consider is 

RJ = ( ' • (30) 
1 o 

In the case of and one has to be careful about the precise definition of 
the specific heat and the reduced temperature. We have checked that, due to a 
cancellation, our definitions of Rg and R^ indeed coincide with those used in the 
literature. 



7.2 Numerical results 

A straight forward way to obtain numerical estimates of amplitude ratios is to fit the 
numerical data for the various observables using ansatze like eqs. (|T1I23II27|29I) . Then 
the amplitude ratios are simply computed from the amplitudes that are obtained 
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from these fits. Here instead we follow a strategy that had already been employed 
in ref. [23] where universal amplitude ratios for the three-dimensional Ising univer- 
sality class had been computed. In order to eliminate the dependence of the result 
on the critical exponents, we consider ratios at a finite reduced temperature (124j) . 
As a first example let us consider 

R T = Yim^ 2nd (-t)T(t) , (t>0) . (31) 

We have computed this product for the reduced temperatures t = j3i ow — j3 c , where 
f3i ow are those values of (3 > j3 c for which simulations have been performed. The 
values for T are taken from tables El [3] and HI The tables [51 [6] and [7| contain no 
exact matches for —t = /3 c — j3 low . Therefore we computed £ 2 ™2 at — t, i.e. at 
fthigh = 2j3 c — Plow, by interpolating the results given in tables EJ [6] and El For this 
purpose, we took two values /3hi a h,i and /3hi g h,2, where we have simulated, such that 
Phigh,i < Phigh < (3hi 9 h,2 and (3high,2-f3high,i is minimal. Then we interpolate linearly 
in fzndiP) = &nd(P)(Pc ~ P) - 6717 to get &nd at f3 high . Our results for &nd{-t)T(t) 
are shown in figure [2j The statistical error of the product is completely dominated 
by the error of T. Unfortunately, the error rapidly increases as we approach the 
critical temperature. We have checked that the error of ^2nd(— t)Y(t) induced by 
the uncertainty of j3 c is negligible. 

Based on RG theory (see e.g. ref. [23]), we expect that the product behaves as 

Ud(-t)T(t) = R r + ht e + Cl t + d x f + b 2 t 29 + ... . (32) 

The leading non-analytic correction is characterized by the exponent 9 = vu). As 
numerical values we take v = 0.6717(1) and uo = 0.785(20) given in ref. [ID] . This 
value of io is consistent with e.g. the result of ref. [25] u = 0.789(11) from the 
perturbative expansion in d = 3 and uo = 0.802(18) from the e-expansion. There is 
little information on 9' in the literature. We assume the value 9' ~ 29 as obtained 
in ref. [26]. Since 29 ~ 1, it is difficult to disentangle the three terms Cit, d\t 9 and 
b2t 29 . Therefore we have used the ansatz 

6nd(-t)T(t) =Ry + bt 9 + Ct (33) 

for our fits. Since A = 2.1 and 2.2 are close to A* we expect that b is small for these 
values of A. We have fitted our results for the 4 model at A = 2.1 and the XY 
model with the ansatz (1331 . In both cases, all available data points were included 
into the fit. In the case of A = 2.2 we have just averaged the available results, since 
we have only results for small values of t. The results of these fits are given in 
table El In all cases we get an x/d.o.f. close to 1. For A = 2.1, the coefficient b of 
corrections oc t e is within error-bars consistent with zero, confirming that A = 2.1 is 
close to A*. Within two standard deviations the result for R r is consistent among 
the three different models. As our final result we quote 

R r = 0.411(2) , (34) 
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Figure 2: The product ^ndi.— t)Y(t) for the </> 4 model at A = 2.1, A = 2.2 and the 
XY model is plotted as a function of t = (3—(3 c , where we have taken f3 c = 0.5091503, 
0.5083355 and 0.4541655 for the <p A model at A = 2.1, A = 2.2 and the XY model, 
respectively. The dotted lines should only guide the eye. 



Table 8: Fits of the product ^2nd(—t)T(t) using the ansatz ( f33)) in the case of the 
<f) 4 model at A = 2.1 and the XY model. In the case of the <j) 4 model at A = 2.2 we 
iave just averaged over all data that are available. 



model 


Ry 


b 


c 


x/d.o.f. 


A = 2.1 


0.4118(8) 


0.001(10) 


-0.28(2) 


1.12 


A = 2.2 


0.4092(11) 






0.67 


XY 


0.4097(18) 


-0.186(28) 


-0.43(8) 


0.22 
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Figure 3: We plot the product g 2 ^"ipsft) ^ or ^ e ^ model at A = 2.1, A = 2.2 and 
the XY model as a function of t = p — p c , where we have taken (3 C = 0.5091503, 
0.5083355 and 0.4541655 for the 4 model at A = 2.1, A = 2.2 and the XY model, 
respectively. The dotted lines should only guide the eye. 

where the error-bar is chosen such that the result for each of the models is covered. 

Next we have computed the amplitude ratio Rb- Also here we have computed 
the ratio of observables at finite reduced temperatures: 

R B = \im 2 * ( ~3^ , (t>0) • (35) 

To this end we have followed a similar approach as for Ry. We took the values 
^(fliow) given in tables [2j [3] and |4] for the 4 model at A = 2.1 and A = 2.2 
and the XY model, respectively Then we computed ^ndiflhigh) and xiPhigh) at 
Phigh = 2/3 c — /3i ow by interpolation of the results given in tables El [6] and for the 
(/> 4 model at A = 2.1 and A = 2.2 and the XY model, respectively. The interpolation 
is done in an analogous way as discussed above. Our results are plotted in figure [31 

In principle, a similar ansatz as for Ry could be used here. However it turns 
out that the analytic background of the magnetic susceptibility has to be taken into 
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Table 9: Fits of the <p A model at A = 2.1 and the XY model using the ansatz / f36j) . 
In the case of the XY model all available data are used, while in the case of the 
cj) 4 model all data up to P = 0.55 are used; i.e. j3 = 0.58 is skipped. In these fits 
j3 c = 0.5091503 and (3 C = 0.4541655 is used as input for the 4 model at A = 2.1 
and the XY model, respectively. 



model 


Rb 


b 


c 


d 


X 2 /d.o.f. 


A = 2.1 


2.828(5) 


0.14(20) 


7.2(2.1) 


14.3(5) 


1.44 


XY 


2.801(15) 


3.2(5) 


-3.1(4.7) 


46.4(7.6) 


2.21 



Table 10: Fits of the 4 model at A = 2.1 and A = 2.2 using the ansatz ( f37p . In 
the case of A = 2.1 we have taken into account all data up to (3 = 0.513, while for 
A = 2.2 all available data have been used. 



model 


Rb 


c 


X 2 /d.o.f. 


A = 2.1 
A = 2.2 


2.8329(26) 
2.8276(29) 


9.7(9) 
13.4(1.0) 


0.24 
1.13 



account to get acceptable fits for a large range of /3-values. Hence we have used the 
ansatz 

2 *(-t) _ o . u+6 



( , =R B + bf + ct + dr . (36) 

The results of these fits are summarized in table M To check the uncertainty due to 
the error of (3 C , we have repeated the fits for shifted values of (3 C (both for computing 
the combination g 2 as wen as m ^ ne ansatz). We get Rb = 2.815(3), 

2.797(16) using (5 C = 0.5091509 and (3 C = 0.4541665 for the 4 model and the XY 
model, respectively. 

We also performed fits for A = 2.1 and 2.2 taking only analytic corrections into 
account: 

e 2 ^%m =RB + ct ■ ^ 

6nd(-^) 3 S(t) 

The results of these fits are given in table [TDl Also here we have repeated the 
fits with shifted values of (3 C . We get R B = 2.83040(26) and 2.82494(29) using 
(3 C = 0.5091509 and f3 c = 0.5083362 for A = 2.1 and A = 2.2, respectively. Based on 
the fits with the ansatz (|37|) we arrive at our final result: 

R B = 2.83(1) . (38) 

The error which is quoted covers the statistical error, the error due to the uncer- 
tainty of p c and the error due to residual leading corrections. The latter is estimated 
by the difference of the results for A = 2.1 and A = 2.2. 



17 



Finally we have computed R£ and . To this end we have analysed our 
very accurate data for the energy density. In the neighbourhood of the critical 
temperature, the energy density behaves as 

E{t) = e ns + c ns t + a±t l - a + ... (39) 

with a±(l — a) = A±. e ns and c ns are the non-singular contributions to the energy 
density and the specific heat at the transition temperature. Their numerical values 
had been accurately determined for the 4 model at A = 2.1 and 2.2 in ref. [6] using 
finite size scaling at the transition temperature. These results are given by [6]: 

e ns = 0.913213(5) + 20 x (& - 0.5091503) + 5 x 10~ 7 x (1/a + 1/0.0151) (40) 

for A = 2.1 and 

e ns = 0.913585(5) + 20 x (f3 c - 0.5083355) + 6 x 10~ 7 x (1/a + 1/0.0151) (41) 
for A = 2.2. The results for the non-singular part of the specific heat are [6]: 

c ns = 157.9(5) + 147000 x (/3 C - 0.5091503) - 2.1 x (1/a + 1/0.0151) (42) 
for A = 2.1 and 

c ns = 155.6(4) + 121000 x (/3 C - 0.5083355) - 2.1 x (1/a + 1/0.0151) (43) 

for A = 2.2. Corresponding results for the XY model are not available. 

Using these results we have computed the singular part of the energy density as 

e s ((3) = E(P) - e ns - c ns (P ~ &) , (44) 

where the numerical values for E(/3) are taken from tables [5] and [6] for the high 
temperature phase and from tables [2] and [3] for the low temperature phase. First 
we have computed 

R+ = ]im& nd {t)[-a{l-a)e s {t){-t)} 1/3 , (t < 0) , (45) 

for the high temperature phase. Note that the combined critical exponent of the 
right hand side vanishes due to the hyperscaling relation dv = 2 — a. Our results 
for this combination are given in figure HI In improved models, the corrections that 
are clearly visible, are to leading order due to terms oc t 2 and oc t 2 ~ a missing in 
eq. (139]) . Our data do not allow to disentangle these two terms. Therefore we just 
linearly extrapolated our results for £,2nd(t) [—Oi(l — a)e s (t)(-t)] 1/3 tot = 0. We get 
from the extrapolation of the 5 largest values of /3 the results Rt = 0.35616(11) and 
i?£ = 0.35626(6) for A = 2.1 and A = 2.2, respectively. The error of R^ is actually 
dominated by the error induced by the uncertainty of our input parameters, e ns , 
c ns , P c and a. In order to estimate this error, we have repeated the whole procedure 
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0.3565 
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* 0.3555 - 
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0.3545 
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Figure 4: In this figure we show R^if) = ^ndif) [— a(l — a)e s (t)(— t)] 1 ^ 3 for the 4 
model at A = 2.1 at A = 2.2. This product is plotted as a function of —t = [3 c — (3, 
where we have taken f3 c = 0.5091503 and 0.5083355 for the 4 model at A = 2.1 
and A = 2.2, respectively. e ns and c ns are set to their central values, v = 0.6717 is 
assumed. 
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for shifted values of these input parameters. In particular, we have replaced e ns by 
e ns + error and similarly for the other input parameters. The errors obtained this 
way are very similar for the two values of A. The largest contribution to the error 
originates from the uncertainty in c ns followed by (3 C , a and e ns . Note that the 
relatively small error due to the uncertainty of a is due to a cancellation of the 
variation of the a that appears explicitly in the definition of Rt and that due to 
the dependence of c ns (I42|43l) on a. Adding up all these errors we arrive at our final 
estimate 

R+ = 0.3562(10) . (46) 

Finally we computed 

R^ = lim ^ [-a(l - a)e.(t) (-t)} 1/3 , (t > 0) . (47) 

Note that in the low temperature phase e s (t) is negative, hence the product con- 
tained in [] is positive. Our numerical results for the 4 model at A = 2.1 and 

2.2 are shown in figure [51 Unfortunately the statistical error blows up for small 
reduced temperatures. This is mainly due to the fact that the relative error of 
the helicity modulus T rapidly increases as the critical temperature is approached. 
In the case of A = 2.1, the value of R^ stays essentially constant over the range 
0.01 < (3 — p c < 0.03. Therefore we regard the apparent increase of the value at 
small (3 — (3 C as a statistical accident. This is confirmed by the fact that for A = 2.2 
no such trend can be seen. Motivated by this consideration, we take our final result 
for RJ as the average over all data with j3 — j3 c < 0.005. We get RJ = 0.852(2) 
and 0.848(3) for A = 2.1 and A = 2.2, respectively. As our final result we take the 
average over the two values of A: 

R^ = 0.850(5) . (48) 

The error-bar is chosen such that the results of both values of A, including their 
individual error-bars, are covered. We have also checked the possible error due to 
the uncertainty of the input parameters a, j3 c , e ns and c ns . Here, in contrast to Rg, 
these errors can be ignored. Actually the result for RJ is virtually independent on 
a. Using the experimental value a = —0.0127 [3l H] instead of our theoretical one, 
the result for R^ changes only in the fourth digit. This is due to the fact that the 
variation of a that appears explicitly in the definition of Rt essentially cancels with 
that of c ns (H2E3D. 

7.3 Comparison with other theoretical and experimental 
results 

In table [11] we have summarized results from the literature that where obtained 
by different theoretical methods and from experiments on the A-transition of 4 He. 
Essentially we followed table 22 of ref. |23j . 
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Figure 5: In this figure we give (t) = T(t) 1 [— a(l — a)e s (t) (— £)] ' for the 4 
model at A = 2.1 at A = 2.2. It is plotted as a function of t — (3 — (3 C , where we 
have taken f3 c = 0.5091503 and 0.5083355 for the 4 model at A = 2.1 and A = 2.2, 
respectively. e ns and c ns are set to their central values, v = 0.6717 is assumed. 
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Table 11: We summarize results obtained by using different theoretical methods and 
experiments on the X-transition of 4 He. MC denotes Monte Carlo simulations of 
lattice models, HT the high temperature series expansion of lattice models, IHT-PR 
the high temperature series expansion of improved lattice models combined with 
the parametric representation of the equation of state, d=3 exp the Geld theoretic 
expansion in three dimensions fixed and e exp the e-expansion. For convenience, 
we have added our present results in the second column of the table. A detailed 
discussion is given in the text. 





MC 


HT, IHT-PR 


d=3 exp 


e exp 


Experiment 


Rt 


0.3562(10) 


0.355(3) [TJ 


0.3606(20) [Ml [27] 


0.36 [IS] 






0.349(5) [T6] 


0.361(4) [30] 








R- 


0.850(5) 




0.815(10) [3TJ[32] 


1.0(2) [IS \M H] 


0.85(2) [34] 




1.180(36) [16] 










Rt 


0.128(2) 


0.127(6) [H] 


0.123(3) [35] 


0.106 [36J[37J[38] 






0.128(4) [IS] 










Rx 


0.411(2) 






0.27 [33] 


0.39 gsmi [331 




0.293(9) [IB] 






0.33 US [3H] 


o.4i memos] 



For most theoretical results are in reasonable agreement among each other. 
In the case of the field theoretic expansion in three dimensions [27] [28] there is 
a discrepancy with our result that is somewhat larger than the combined error. 
There is no experimental result for this amplitude ratio, since there is no direct 
experimental access to the correlation length of 4 He in the high temperature phase. 
The authors of ref. [32] quote their Monte Carlo result as a function of a: = 
0.3382(14) - 0.717(96)a + 0.87(1. 13)a 2 . We have inserted a = -0.0151 to get the 
value quoted in table [TTJ It differs by about 3 times the combined error from our 
result. 

In the case of R^ our result is in perfect agreement with the experimental number 
given in table IV of ref. [31]. It is interesting to note that also the experimental value 
of Rj shows very little dependence on the value of a that is assumed in the analysis; 
the authors of ref. (33] arrive at R7 = 0.86 and R^ = 0.84 using a = —0.007 and 
a = —0.025, respectively, instead of their preferred value a = —0.016. They do 
not explicitly quote an error for their estimate of R7- In section III.B. they write 
however that {R^) 3 is constant within 5% over the entire range of pressures from 
to 30 bars. In the review [2] the result of ref. [33] is quoted as R7 = 0.85 ± 0.02. 
Concerning theory, there is a discrepancy by about twice the combined error with 
the field theoretic result [31] [32]. The error of the result obtained by the e-expansion 
[29] [331 [2] is quite large. The authors of ref. [16] quote their Monte Carlo result 
as Rj = 1.1580 — 0.696a + 0.97a 2 ± 0.036. The huge difference compared with our 
result is likely due to the difference in the estimates of T as discussed in section 

In the literature our universal amplitude ratio Rb is not discussed. Instead one 
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finds 




aA + C. 



+ 



(49) 



It can be expressed as 




{RtfR B = 0.128(2) 



(50) 



where the error is dominated by the error of R^. Our present result is in good 
agreement with that obtained in ref. [14] using the high temperature expansion of 
improved models in combination with the parametric representation of the equa- 
tion of state. It also agrees with ref. |16| . where a combination of Monte Carlo 
simulations and the parametric representation of the equation of state was used. In 
ref. [16] i?+ is given as a function of a. The value in our table is at a = —0.0151. 
Again there is a discrepancy with the field theoretic result [35] that is slightly larger 
than the combined error. The e-expansion for i?+ has been computed up to 0(e 2 ) 
[36] 137] 138] . The number given in the table is obtained by simply setting e = 1. 
There is no experimental result for since there is no experimental access in 4 He 
to the analogue of the magnetisation in spin models. 

In the case of Rr we see the largest differences among the results obtained by 
using theoretical methods. The Monte Carlo result, quoted in eq. (97) of ref. [16] , 
differs by more then 10 times the combined error from ours. This huge difference 
can be traced back to the discrepancy in T at given values of (3 as discussed in 
section [5] Hohenberg et al. [33] have computed Ry by using the e-expansion to 
O(e). Bervillier [29] extended the calculation up to 0(e 2 ). Okabe and Ideura [39] 
corrected the calculation of Bervillier, which does however not change the numerical 
value, and computed the 1/iV-expansion resulting in R-? = 0.14 for N = 2. 

There are no direct experimental results for the correlation length in the high 
temperature phase of 4 He. Instead one might use the data for the specific heat in 
combination with the theoretical result for R^ to arrive at the amplitude of the 
correlation length. This way, using R^ = 0.36 and the experimental results of refs. 
[USE], Hohenberg et al. [33] arrive at i?x = 0.39. Bervillier noted, see section 
III. A. of [21], that there is an error in the experimental value of the amplitude of 
the transversal correlation length used by Hohenberg et al. [33]. He arrives at 
the corrected value Rr = 0.41. It would certainly be worth while to redo this 
calculation using most recent experimental data; e.g. those of ref. [3] H] for the 
specific heat and our estimate of R^. 

8 Summary and Conclusions 

In this paper we have computed universal amplitude ratios for the three dimensional 
XY universality class. These results are based on Monte Carlo simulations of the 
three dimensional XY model and the 4 -model at A = 2.1 and A = 2.2. Note that 
these values of A are close to A* = 2.15(5) [ID] , where leading corrections to scaling 
vanish. We performed simulations in the low and the high temperature phase of 
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these models. Extracting results for the thermodynamic limit, one has to take into 
account the effect of the Goldstone mode in the low temperature phase. For a 
discussion see section HI 

Our results for universal amplitude ratios for the three-dimensional XY univer- 
sality class are throughout more precise than previous theoretical results. They 
are obtained in a rather direct way, making hidden systematic errors unlikely. 
When available, our results agree nicely with experimental ones obtained for the A- 
transition of 4 He, giving further confirmation to the fact that this transition shares 
the three-dimensional XY universality class. 

The numerical results obtained here for the correlation length ^2nd, the helicity 
modulus T, the energy density E and the specific heat set the stage also for the 
study of the specific heat or the Casimir force in confined geometries using improved 
models. 
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